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[For the first time, the Fermi-LAT measured the angular power spectrum (APS) of anisotropics in the diffuse gamma-ray 
■background. The data is found to be broadly compatible with a model with contributions from the point sources in 
.the 1-year catalog, the galactic diffuse background, and the cxtragalactic isotropic emission; however deviations are 
present at both large and small angular scales. In this study, we complement the model with a contribution from Dark 
^Matter (DM) whose distribution is modeled exploiting the results of the most recent A-body simulations, considering 

Othe contribution of cxtragalactic halos and subhalos (from Millenium-II) and of galactic substructures (from Aquarius) . 
^Vith the use of the Fermi Science Tools, these simulations serve as templates to produce mock gamma-ray count maps 
for DM gamma-ray emission, both in the case of an annihilating and a decaying DM candidate. The APS will then be 
.computed and compared with the Fermi-LAT results to derive constraints on the DM particle physics properties. The 
Onpossible systematic due to an imperfect model of the galactic foreground is also studied and taken into account properly. 
Q ^he present paper reports on the status of the project. 
■ 
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1. Introduction 

The Isotropic Gamma-Ray Background (IGRB) can be 
defined as the radiation that remains after the resolved 
sources (both point-like and extended) and the galactic 
foreground (produced from the interaction of the cosmic 
jays with the interstellar medium) are subtracted from the 
total gamma-ray emission. The most recent measurement 
of the IGRB energy spectrum was done by the Fermi-LAT 
telescope and presented in Refs. [H, [l^L Contrary to what 
[was found by the EGRET telescope]! Q, the IGRB ap- 
■pears to be perfectly compatible with a power-law with a 
slope of —2.41. at least up to 10 GeV. 

The contribution of unresolved sources, both cxtra- 
galactic (e.g., blazars 0, H, 0, H, S jiol - E| , star- forming 
galaxies [y, [l^ and radio galaxies |l4l Il3| ) and galactic 
sources (like milli-second pulsars [15|) have been consid- 
ered. The contribution of each class is estimated from 
population studies of the detected objects (see Ref. [lH 
for an example in the case of blazars) and it turns out that 
all the components considered up to now may not be able 
to account for the whole IGRB [l^. Thus there may be 
room for additional contributions like, e.g., that of Dark 
Matter (DM) [B, El, i, H H ■ Actually, the authors of 
Ref. have already used the Fermi-LAT measurement of 
the IGRB intensity to constrain the properties of the DM 
particle, finding that, in their most optimistic scenario. 



values larger than 10~^^cm'^s^^ for the annihilation cross 
section ((Ta„„w) can be excluded. 

Complementary to the energy spectrum of the IGRB, 
we can also the recent analysis of the angular power spec- 
trum (APS) of anisotropics in the diffuse gamma-ray back- 
ground [m, [u 



231 to study the nature of the IGRB. The 



analysis was conducted by the Fermi-LAT collaboration 
and resulted in the detection of some angular power above 
multipoles oi i > 155, with a significance larger than Scr, 
in the each of the energy bins between 1 GeV and 10 GeV. 
More detailed information on this measurement will be 
provided in Sec. [5J We simply note here that, even if such 
a detection seems to be compatible to the signal predicted 
for a population of unclustered, unresolved blazars, it can 
still be used to put some useful constraints on the nature of 
the DM particle. The goal of the present work is to derive 
those constraints using state-of-art numerical simuations. 

The project is divided in two parts: in the first one 
we will make use of the most recents results from A-body 
simulations to derive all-sky maps of gamma-ray emission 
from DM annihilations and decay. We will consider the 
contributions of i) cxtragalactic halos and subhalos (based 
on the Millennium II A-body simulation [l^l), ii) the emis- 
sion from the smooth halo of the Milky Way (modeled as in 
Ref. I25I a nd in) its subhalos (from the results of Aquar- 
ius [26|, |27|). Moreover, we will also consider the emission 
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from halos and subhalos below the mass resolution of the 
two simulations mentioned above, down to the minimal 
halo mass M„i„. The all-sky maps will then be used to 
compute the APS of anisotropies in the gamma-ray emis- 
sion from DM annihilation and decaj^ We plan to present 



the results of this first part in Ref. l28| , where we will also 
study the dependence of the APS from the assumptions 
made while building the maps of DM-induced gamma-ray 
emission. Refer also to Sees. [3] and H] for more information 
on this first part. 

The second part of the project will deal with the com- 
parison with Fermi-LAT APS data. The simulation maps 
produced in Ref. |28[ will be used as templates for the 



DM contribution to the IGRB. We will use the most recent 
Fermi Instrument Response Functions (IRFs) to estimate 
how experimental issues may affect the APS. It will be 
particularly important to take correctly into account the 
experimental Point Spread Function (PSF) and the Fermi- 
LAT exposure. We will finally use the APS measurement 
to derive constraints on the nature of DM taking into ac- 
count also the constraints from the new measurement of 
the IGRB energy spectrum (see also Sec. [5]). 

We will conclude this small introduction reminding that 
this note should be considered just as a status report, since 
none of the two parts of the project is complete up to now. 
We will present some preliminary results in Sees. [3] and HI 
but we will mainly focus here on the methodology, empha- 
sizing that the final results will be extensively presented 
in the near future 



28|. 



2. Fermi-LAT measurement of the anisotropies in 
the diffuse gamma-ray background 

For the analysis in Ref. the first 22 months of 

Fermi-LAT data have been analyzed, dividing the energy 
range between 1 GeV and 50 GeV in 4 energy bins. The 
point sources in the first 1 year catalog |29| have been 
masked, as well as the emission within a band of 30 de- 
grees above and below the galactic plane. This was done 
to cover the regions in the sky where the emission is dom- 
inated by resolved sources and by the galactic foregound, 
and to restrict the analysis only to where the IGRB is 
a significant component. We note here that the analy- 
sis performed in Ref. [2^ to measure the IGRB intensity 
spectrum is different since templates of emission are imple- 
mented to subtract known components in the gamma-ray 
emission, instead than masking portions of the sky. Thus, 
strictly speaking, the data used in Ref. [l^ and analyzed 
to derive the APS are not the same used for the IGRB 
energy spectrum in Ref. Moreover, even at high lat- 

itudes (|/| > 30°) the galactic emission is still important 
and non-negligeble. We therefore expect that some level of 
contamination in the APS may be present from this kind 
of background. These contaminations are supposed to be 
located primarly at low angular multipoles £ (large an- 
gular scales) and so only multipoles larger than 155 have 



been considered in Ref. [23j. On the other hand, mul- 
tipoles larger than 504 are also discarded since, at these 
small angular scales, the signal is strongly damped by the 
experimental PSF. 

Two definitions of APS are used: i) the intensity APS 
as in Eqs. [T] and [2] below, for which the intensity gamma- 
ray maps arc decomposed directly in spherical harmonics 



= I /(vi/)y,;„(vi/)df], 



(1) 



(2) 



and a) the fluctuaction APS, that can be derived by 
the intensity one, dividing for the average intensity squared. 

As a consequence, the fluctuation APS will not depend 
on the energy of the emission, if the distribution of the 
sources is energy-indenpendent. 

The Fermi-LAT reported detection of angular power 
in all the 4 energy bins considered, with a significance of 
larger than 3(t in the energy bins from 1 GeV to 10 GeV. 
The data have been compared with the power spectrum 
of a source model made of i) the point-like sources in Ref. 
(jgj . a) a model for the galactic foreground and Hi) an 
isotropic component at the level of the IGRB in Ref. [20| . 
In the same regions outside the mask defined above, the 
general features in the APS of the model are similar to 
those in the data, but the model APS does not accurately 
reproduce the data APS in all energy bins on small or large 
angular scales. Furthermore, the model angular power at 
155 > £ > 504 is consistently below that measured in the 
data. 

This seems to point to the possibility of having a pop- 
ulation of unresolved sources contributing to the APS, in 
order to explain the data. The fact that the intensity 
APS is compatible with being independent of multipole, 
together with the way the intensity and fiuctuation APS 
changes in the 4 energy bins suggests that one or more 
populations of unresolved, unclustered classes of sources 
may be responsible for the missing power. 

If this hypothesis is true, then the data will provide us 
some useful constraints on the DM particle: the normaliza- 
tion and shape of the APS from DM annihilation has been 
studied extensively in the last years, focusing on the case of 



extragalactic halos and subhalos Il8l |9|. | 3 d. l31l. |32| , on the 
ease of a galactic component 3J, |33|, |lO[ or both [19|, |35 |. 
Comparing model predictions with the Fermi-LAT APS 
data allows us to put constraints on how important are 
DM-induced gamma-rays in the IGRB and, consequently, 
draw exclusion lines on the annihilation cross section, or 
other quantities more difficult to constrain like the mini- 
mal halo mass Af,„i„. 

3. Extragalactic halos and subhalos 

This section and the next one will be devoted to de- 
scribe the methodology that will be used to compute the 



2 



maps of DM-induced radiation. To model the extragalac- 
tic emission, we use the catalogs of halos and subhalos of 
the Millennium-II A^-body simulation; the simulation box 
is a cube with a size of 100 Mpc//i and it contains DM 
halos and subhalos down to a mass resolution of M^^ = 
6.89 X 10^ Mq/H Halo catalogs are available at mul- 
tiple snapshots with different redshifts up to z = 127. We 
analyze them in a very similar way to what has been done 
in Ref. [3^ : we will only consider objects with more than 
100 particles and derive the properties of each halo (DM 
profile, luminosity and concentration) from the maximal 
circular velocity Knax and radius r^^^ (corrected from spu- 
rious effects related to numerical resolution, as done in 
Ref. [321) and assuming a Navarro-Frenk- White (NFW) 
profile |36|. Since we want to probe a volume which is 
much larger than the Millennium-II simulation box, we 
implement the same technique described in Ref. [32| to 
construct sky-maps of the extragalactic signal, dividing 
the past light-cone in concentric shells (each of them at a 
particular redshift) and then filling them up with identi- 
cal copies of the Millennium-II simulation boxes at that 
particular redshift (see Fig. 9 in Ref. [13 )• Each cube 
is randomly rotated and translated in order to avoid any 
boundary effect. In this way one can compute the lumi- 
nosity from a particular direction in the sky simply by 
accounting for all the halos encountered in the direction of 
This was done for all halos above Millennium-II mass 
resolution and you may see the results in Fig. [T] up to a 
z = 2.6. 

We also want to consider the contribution of halos and 
subhalos below the mass resolution of Millennium-II. We 
treated main halos and subhalos separetely and in two dif- 
ferent ways. We will start by describing how we account 
for main halos with a mass smaller than M^^^. Actually, 
the authors of Ref. [s^] already considered this contri- 
bution: fitting the Millennium-II halos, they derived the 
cumulative luminosity function F{M) = ^L/AfAlogM, 
providing the total luminosity of main halos with a mass 
between M and M -\-dM . Then, they assumed that F{M) 
can be extrapolated below M^os in order to compute the 
gamma-ray flux predicted from main halo between M„i„ 
and M„5. This quantity was then used to boost up the lu- 
minosity of halos with a mass between 1.4 x IO^A/q/Zi and 
6.89 X IO^Mq/Zi. This was done under the assumption that 
the distribution of main halos below A/^^ follows exactly 
that of those between 1.4 x IG^Afo/Zi and 6.89 x lO^Afg/Zi, 
such an approach is by the fact that the two-point correla- 
tion function (which is an indicator of clustering) is found 
to reach a constant value when approaching Af^^ (see Fig. 
10 of Ref. HI). 

We want to test this assumption by implementing a 
different way of including the contribution of halos below 
Af,„5 . We will then look for differences in the APS to check 
if the approach used to cover the range between Af„i„ and 
M^^^ affects in any sense the shape of the APS. We there- 
fore proceeded as follows: 



• divide the range between M^-,^ and Af„B in mass 
decades (Af^, Af,;+i). For each mass decade, we con- 
sidered all the halos in the Millennium-II simula- 
tion box with a mass between 1.4 x lO^Afo/Zi and 
6.89 X 10^A^q//i and assign to each of them a new 
mass between Af,; and Mij^i, assuming masses fol- 
low a probability distribution equal to the halo mass 
function dn/dM. 

• we also assign a luminosity to each of these halos. 
From these two quantities, we can completely derive 
the halo profile. In this way, we end up by having 
a box of simulated objects with masses between Mi 
and Afi+i distributed as those between in the range 
1.4 X 10^ - 6.89 X IO^Mq/H. 

• however, in the mass decades we are considering, the 
halo mass function will predict more halos than we 
have actually have in the simulated box. So, we stack 
together different copies of the box obtained in the 
previous point until we reach the correct number of 
halos. The different copies are randomly rotated one 
with respect to the other. 

• the box that comes out from the stacking is used 
to fill up the region up to a maximal distance of 
i?„ax, defined as the distance above which all halos 
in a given mass decade become point-like. Maps are 
produced from this distribution of objects till i?„ax- 

• for the regime beyond i?„ax, we do not consider any 
stacking, we simply take the boxes coming from point 
2 of this list and create the sky- map from those. Mul- 
tiple copies of the map are then stacked together till 
the dcsidered flux is reached. Since we are now in the 
regime where halos are point-like, the procedure is 
similar to what has been done below R^^^ but with- 
out having to consider the 3D distribution of halos. 

This procedure, in a sense, moves in the opposite di- 
rection to what has been done in Ref. [IH since we start 
by assuming that halos below Af^^^ have the same distri- 
bution of those above, but then we pass through a lot of 
independent rotations so that, at the end, the final maps 
are expected to be more isotropic than those in Ref. ^32i] . 
Once we compute the APS from the two approaches we 
will be able to see if differences show up in the shape of 
the APS. 

In Fig. [5] (left panel) the black crosses indicate the 
IGRB energy spectrum taken from Ref. The red 

points show the amount of flux produced in halos and sub- 
halos above M^^^ until z = 2.6 (those depicted in Fig. [ij 
in the case of annihilating DM (filled points) and decaying 
DM (empty points). If we also consider the emission be- 
low Af,„5 (implemented as described in the list above) the 
emission increases to be that one of the green points. For 
annihilating DM the increase is approximately of an order 
of magnitude, while for decaying DM it is completely neg- 
ligible (red and green empty points practically overlap). 
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Figure 1: All-sky maps of gamma-ray emission from DM annihilation (left panel) or decay (right panel) in halos and subhalos above the 
mass resolution of Millennium-II, up to 2: = 2.6. Emission comes from hadronization of b quarks and is computed at 10 GeV. For the case of 
annihilating DM, the mass is 200 GcV and the annihilation cross section 3 X lO^^^'cm^s^^ , while, for the decaying case, the mass is 2 TeV 
and the decay lifetime is r = 2 X lO^^s. 



The total flux is still a factor of 100 (50) for annihilation 
(decay) below the IGRB, at least at 10 GeV. 

In order to complete the description of the cxtragalac- 
tic component, the only remaining ingredient is the emis- 
sion from subhalos below M^^. Those will be described 
following the prescription derived in Ref. (STj and gener- 
alized in Ref. [SSj . In the first reference the subhalos of the 
Via Lactea II A^-body simulation are analyzed to compute 
the probability distribution P{p, r) of having a DM den- 
sity between p and p + dp a.t a distance r from the center 
of the MW halo. P{p, r) has two distinct components: i) a 
parabolic regime due to the smooth halo and ii) a power- 
law that extends to larger densities due to substructures 
(see Fig. 1 of Ref. (s^])- The authors, then used P{p,r) 
to coinpute the boost factor due to substructures, while 
Ref. [38[ extended the prescription to halos larger and 
smaller than the MW halos, having in mind the case of 
galaxy clusters and dwarf Spheroidal galaxies respectively. 
We will use it to boost up the emission of the halos below 
M,35, , accounting in this way for the emission of unresolved 
subhalos. 

4. The Milky Way halo and its subhalos 

We continue now with the implementation of the gamma- 
ray emission related to our own galaxy: we use the distance 
of 700 kpc (approximately 3 times the virial radius of the 
Milky Way (MW) halo) as the limit between the extra- 
galactic and the galactic regime. 

The smooth halo of the MW is modeled with a NFW 
profile with parameters taken from Ref. [25j , a model that 
is consistent with the available observation data on the 
MW. The total flux of the MW smooth halo is plotted in 
the left panel of Fig. [51 It represents the largest among 
the different DM components plotted there. However we 
should note that at the moment of computing the APS, 
we will mask the galactic plane, as it has been done by 
the Fcrmi-LAT collaboration in Ref. This will have 



the effect of drastically decreasing the emission associated 
with the smooth halo. 

To include the galactic subhalos we will use the Aquar- 
ius A^-body simulation . The same procedure described 
in the previous section will be used here for the Aquarius 
subhalos, without the need to considering replicas of the 
simulation box since we are interested only in the region 
within 700 kpc. The emission associated with galactic sub- 
halos is indicated as yellow points in Fig. [21 We have not 
yet implemented the contribution of subhalos below the 
mass resolution of Aquarius M/^^ = 1.712 x 10"^ Af©. This 
will substantially increase the emission associated with 
galactic subhalos, which is only subdominant if we only 
consider objects above Af'^^. 

5. Deriving constraints 

The maps that will be obtained from the procedure 
sketched in the previous sections will then be used to com- 
pute the APS and to compare the results with the Fermi- 
LAT data. Some results can be seen in the central and 
right panels of Fig. [51 where the total fluctuation APS is 
plotted (black line). Also the contribution of the differ- 
ent components is present, multiplied by the square of the 
average flux of each component with respect to the to- 
tal average emission, so that the sum of the colored lines 
gives the total (black line) also visually. Let us stress that 
these are very preliminary results: not all the components 
have been included (the contribution of unresolved extra- 
galactic and galactic subhalos is missing) and we are not 
considering any mask. 

We can, however, already see some interesting trends: 

• the galactic component dominates the case of annihi- 
lating DM. This may change after we will include the 
contribution of galactic subhalos less massive than 

• the extragalactic component below Af^^^ (both for 
annihilating and decaying DM) does not indicate any 
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Figure 2: Left panel: Black crosses indicate the energy spectrum of the IGRB taken from Ref. [20II . Filled (empty) points refer to the case of 
annihilating (decaying) DM. Red points indicate the flux coming from DM halos and subhalos above Mr^s, while green points only consider 
the extrapolation down to M^i„ for main halos. Blue points show the emission from the smooth DM halo of the Milky Way, while the yellow 
ones accounts for the emission from all the subhalos in the Aquarius simulation. The grey lines refer to the annihilation and decay energy 
spectrum (from the hadronization of b quarks) normalized to blue points. Central and right panel: APS of anisotropics in the gamma-ray 
emission from DM annihilation (central panel) and decay (right panel). The total fluctuation APS is plotted in black, while the colored lines 
correspond to the different components, multiplied by the square of their average flux with respect to the total average flux. The red line 
indicates the contribution of extragalactic halos and subhalos above the mass resolution of Millennium II. The green line refers to main halos 
below Mres and down to M„in, while the blue one accounts for the subhalos in the Aquarius simulation. 



intrinsic clustering, as expected. 

• the contributions of Aquarius subhalos (blue lines) 
and of Millennium-II halos and subhalos (red lines) 
decrease with multipole, a consequence of the fact 
that the APS is sensitive to the inner structures of 
the halos. 

In Fig. [2] we have also plotted the Fermi-LAT measure- 
ment of the fluctuation APS in the energy bin between 1 
GeV and 2 GeV. Some other ingredients are still needed 
in order to be able to perform a proper comparison be- 
tween the DM APS and the data. This will be done in 
the future, taking into account experimental features like 
the effect of the PSF, or a possible residual contamination 
from the Galactic foreground. We plan to include the IRF 
of the telescope using the Fermi Science Tools and treating 
the DM component in the same way as one of the three 
different components of the source model described in Sec. 

m 



APS by the Fermi-LAT collaboration in order to be able 
to put constraints on the nature of the DM particle. As 
the project is still on-going, we only presented here the 
methodology and some preliminary results. We refer in- 
terested readers to keep an eye on the arXiv webpagc for 
the publication of the final results [l^ . 
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6. Conclusions 

Here we briefly summarized a project which is still in 
progress. The goal is to compute realistic and complete 
all-sky maps of gamma-ray emission from DM annihilation 
and decay, taking into account the contribution of both ex- 
tragalactic and galactic halos and subhalos. For the most 
massive objects we will refer directly to the results of two 
of the most recent iV-body simulations (Millcnnium-II and 
Aquarius), but we will also account for smaller halos, down 
to their minimal mass M„,i„. Once obtained, these maps 
will be used to compute the APS of anisotropics. Our 
goal is to compare it with the recent measurement of the 
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